Spatial Coherence Resonance near Pattern-Forming Instabilities 
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The analogue of temporal coherence resonance for spatial degrees of freedom is reported. Specif- 
ically, we show that spatiotemporal noise is able to optimally extract an intrinsic spatial scale in 
nonlinear media close to (but before) a pattern-forming instability. This effect is observed in a 
model of pattern-forming chemical reaction and in the Swift-Hohenberg model of fluid convection. 
In the latter case, the phenomenon is described analytically via an approximate approach. 
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PACS numbers: 05.40.-a, 05.45.-a 

The ability of noise to induce temporal coherence in 
nonlinear systems is a well documented fact. In stochas- 
tic resonance (SR), for instance, random fluctuations en- 
hance the response to weak periodic driving, as observed 
in many different physical, chemical and biological sce- 
narios 1] . From a more intriguing perspective, even sys- 
tems showing no explicit time scale exhibit an enhance- 
ment of temporal coherence due to noise 0, 0. This 
autonomous stochastic resonance, where noise extracts a 
hidden, intrinsic time scale of the system's dynamics, has 
been termed coherence resonance (CR) 0, and has been 
predicted theoretically in a wide variety of models (see Q 
for a review) and observed experimentally in fields as di- 
verse as plasma physics Q , laser dynami cs H , electronics 
[3, atom trapping [j| and neuroscience |ltj|. 

Parallel to the advances in the fields of stochastic and 
coherence resonance, recent years have witnessed an in- 
creasing interest in the dynamics of systems with spa- 
tial degrees of freedom. Much work has been devoted 
to analyze the enhancement of spatiotemporal order by 
noise |llj, |l2( , and the phenomena of standard temporal 
SR and CR have been studied in detail in arrays of dy- 
namical elements |0, Q] , where it has been shown that 
coupling enhances the noise-induced temporal coherence. 
But so far, up to our knowledge, no analogues of CR (or 
SR, for that matter) for spatial degrees of freedom ex- 
ist. Namely, we ask ourselves if noise is able to extract 
and optimize an intrinsic spatial scale from a determinis- 
tically homogeneous extended medium. In this Letter we 
show that such spatial coherence resonance (SCR) exists 
close to pattern- forming instabilities, mimicking one of 
the mechanisms of standard (temporal) CR. 

There arc different mechanisms giving rise to tempo- 
ral CR. In the most extended situation, CR arises as the 
result of a matching between the time scale associated to 
noise-induced escape from a stable steady state and the 
characteristic time scale of the phase-space trajectories 
outside from the immediate vicinity of the fixed point 
(such as what happens in excitable media 0,^3 or close 



to a saddle- node bifurcation in phase models Q). A sec- 
ond, less studied mechanism takes place near the onset 
of dynamical instabilities giving rise to periodic behavior, 
where noise excites precursors of the bifurcation in the 
form of stochastic limit cycles with a well-defined peak 
in the power spectrum, even before the bifurcation really 
happens Il7|. Such noisy precursors exhibit an op- 
timal coherence at an intermediate noise strength, in a 
way typical of coherence resonance [l^. This effect has 
been confirmed experimentally |l9j| . 



Analogously to the temporal case, noisy precursors of 
pattern-forming bifurcations have been observed, in the 
form of noise-induced spatial structures, in electrocon- 
vection J2jj, Rayleigh-Benard convection [2l| and optics 
22], for instance. In this case, the spatial power spectrum 
(now structure function) of the system exhibits a well de- 
fined peak due to noise, even below the bifurcation point. 
In what follows, we show that noise controls the shape 
of this spectral peak in a manner characteristic of coher- 
ence resonance, but with the role of temporal frequency 
being played by spatial wave number. This allows us to 
define a pure spatial analogue of temporal coherence res- 
onance. The generality of the phenomenon is illustrated 
through its occurrence in two different pattern-forming 
systems, namely an activator-inhibitor model of chem- 
ical pattern formation and the Swift-Hohenberg model 
of fluid convection. In the latter case the effect is inter- 
preted analytically. 



We first illustrate the phenomenon of SCR in chem- 
ical pattern formation, in particular in the case of the 
chlorine dioxide-iodine-malonic acid (CDIMA) reaction, 
which is well known to exhibit both temporal and spa- 
tial dynamics |2^| . Experimental and theoretical studies 
of this reaction have concluded that the malonic acid 
serves only as a source of iodide, and thus the system 
can be described by just two variables, in the way of an 
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activator-inhibitor model: |23j, |2, 
ii = V 2 u + a — cu — 4 

v = b i V 2 w + 



(1 + u 2 ) 

U V 



(1 + u 2 ) 



(1) 



where u and v play the role of activator and inhibitor 
fields, corresponding to the concentrations of iodine and 
chlorine dioxide ions, respectively (all variables and pa- 
rameters are in dimensionless units). This model de- 
scribes the photosensitivity of the CDIMA reaction, with 
<f> accounting for the illumination level, which is an ex- 
ternally controlled parameter. Additionally, for large dif- 
ferences in the diffusion coefficients of the two species 
this system presents a Turing instability. In particu- 
lar, the model has an homogenous stable state given by 
Uh = (a — 50)/(5c) and Vh = a(l + u 2 l )/(5uh), which be- 
comes unstable versus spatial perturbations of nonzero 
wavenumber for <f> < cf> c . For a — 16, b = 301, c = 0.6 
and d = 1.07, the critical illumination is <fi c = 2.3 |2J]. 

This chemical system allows for an experimental study 
of stochastic effects, by introducing noise in the illumina- 
tion profile. In this way, for instance, it has been shown 
experimentally that a quenched spatial noise sustains 
Turing patterns in this system well beyond the determin- 
istic bifurcation threshold |2^. We now consider a spa- 
tiotemporally stochastic illumination profile, <fi — > <f> + £, 
where £ (x, t) is a Gaussian noise with zero mean, white 
in time and with spatial autocorrelation given by 

(£{x, t)£(&, 0) = 2 <r 2 C(x - x') S(t - t'). (2) 

In order to examine the spatial response of the system in 
the presence of noise, we evaluate the structure function 
of the activator variable, defined as 



S(k,t) 



(u(k,t) u(-k,t)) 
V 



(3) 



where u(k, t) is the spatial Fourier transform of the ac- 
tivator field, V is the d-dimensional volume of the sys- 
tem, and (• • • ) represents an ensemble average over noise 
realizations. Taking into account the spherical symme- 
try of this system, we compute the spherical average of 
the structure function as s(k) — J n dflk S(k,t), where 

k = k and flk is a hyperspherical shell of radius k. The 

result is shown in Fig.^a) for increasing noise intensities 
|25l |. evaluated for a two-dimensional system and in the 
case of spatially white noise, i.e. C(xi — Xj) — Sij/Ax. 
The average illumination to which the system is sub- 
jected is above threshold, thus for vanishing or very weak 
noises the system is homogeneous and the structure func- 
tion does not show a well-defined peak as k increases. 
However, increasing the noise intensity induces a peak 
in the structure function at k = /c max , which represents 



a spatial noisy precursor of the pattern-forming bifurca- 
tion that occurs at 4> = 4> c : fluctuations are capable to 
anticipate the instability, making the structure manifest 
even before the threshold has been crossed. If fluctua- 
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FIG. 1: (a) Spherical average of the structure function for 
increasing noise intensities. Average illumination is <j> = 2.35. 
(b) AS versus noise intensity for different values of the average 
illumination: <j> = 2.32 (squares), <f> = 2.35 (circles), and (j> = 
2.37 (diamonds). The simulations have been performed on a 
square lattice of 64x64 cells with mesh size Ax = 0.5. 

tions become too large, disorder inevitably comes into 
play, and the peak in the structure function becomes less 
well defined [see triangles in Fig. HJa)]. It is then clear 
that an optimal noise intensity exists at which the struc- 
ture function peak is best resolved from the background 
fluctuations. In order to quantify the efficiency of the 
noise, an adequate measure must be used. There are 
different ways of quantifying such a response depending 
on the peculiarities of the problem under study. Here, 
we estimate the spatial coherence of the system through 
the quantity AS = s(k max ) — s(0), taking into account 
that s(0) measures the level of fluctuations existing in 
the system. Figure 0Jb) shows how the SND varies with 
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FIG. 2: Spatial profiles of the activator field for increasing 
noise strength and cp — 2.35, with color mapped nonlinearly 
(with a tanh function) and the same contrast in all figures 
for comparison. From left to right, a 2 = 0.001,0.2,0.6. The 
lattice parameters are those of Fig. 

noise intensity a 2 for different illumination levels. It can 
be seen that there is an optimal noise strength a 2 pt for 
which the structure function peak is resolved best. This 
is the signature of a spatial coherence resonance driven 
by the external stochastic forcing, which extracts the in- 
trinsic spatial distribution of the system. As shown in 
Fig. nfb), the closer the system is to the bifurcation point, 
the more pronounced the resonance is, and the position 
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of the maximum cr 2 pt shifts to smaller values. 

Figure |2] shows several patterns produced by model 
||TJ for different values of the noise strength at <j) — 2.35. 
While the system is quasi-homogenous for small noise 
(left plot), it exhibits a rather self-organized configura- 
tion for a 2 ~ cr 2 pt (middle plot), and becomes gradually 
more disordered for large noise strengths (right plot). Al- 
though the effect may not appear very clear in the spatial 
profiles, the relevant fact is the appearance and enhance- 
ment of the peak in the structure function, since this is 
a measurable effect through, for instance, scattering ex- 
periments. 

We may illustrate analytically this phenomenon in the 
case of a simpler pattern-forming model, namely the 
Swift-Hohenberg equation, which describes the onset of 
hydrodynamic convection in Rayleigh-Benard cells |26| : 



d(j)(x, t) 
dt 



= - r< j,- < j? -(V 2 + fc 2 ) 2 <j> + i{x,t), (4) 



where the control parameter r is proportional to the tem- 
perature gradient driving the fluid, and the additive noise 
£(x, t) is again considered Gaussian and white in both 
time and space. As in the CDIMA model, in the absence 
of noise and for r > 0, this system has a homogeneous 
steady state <f> = 0, representing conduction, which upon 
decrease of r becomes unstable at r = r c = versus static 
perturbations of nonzero wavenumber k — ko, which rep- 
resents convection |2(t| . In the presence of a small amount 
of noise, the structure function can be calculated by lin- 
earizing Eq. Q around 6 = 0, and transforming the 
system to Fourier space llj. The result is: 



S st (k) 



(fc 2 - k 2 ) 2 + r 



(5) 



Since we want to operate in regimes where, even though 
the system is dcterministically homogeneous, noise can 
be large, we generalize the previous expression J5J by al- 
lowing the parameters cr and r to depend on noise; we 
will call them a' (a) and r'(cr) in what follows. In order 
to estimate how these renormalized parameters depend 
on the noise intensity, we rederive expression (jSJ) with- 
out neglecting completely nonlinear terms. To do that, 
we estimate the cubic term in Q by applying the Gaus- 
sian approximation [27[ 4)(x,t) 3 ~ 3^(x, t)(4> 2 ) st . This 
renders back the equation linear, so that the structure 
function obeys again expression JSJ with r being replaced 
by r' = r + 3((f> 2 ) st . We can now calculate the value of 



the average squared field 
relation 



from the self-consistency 



(2tt)s 



S s t(k)2Ttk dk 



7T 



k 2 



Since we will be operating close to threshold, r 
be small and r' ~ 3(</> 2 ) S f, which leads from £jj 



(6) 

will 

to 



{4> 2 )st ~ (cr 2 ) 2 / 3 . Therefore, from the definition of r' 
above we finally obtain that this effective parameter also 
scales with noise intensity as r' ~ (cr 2 ) 2 / 3 . In order to 
verify this result, we perform numerical simulations of 
model (@J) in the homogeneous regime for different noise 
intensities, and fit the resulting structure functions with 
expression ||5}, with effective parameters a' and r' which 
are allowed to depend on noise. The results of this anal- 
ysis are shown in Fig. and confirm that cr' scales lin- 
early with noise intensity a (as expected, since nonlinear 
corrections did not affect this parameter in the analysis 
made above), whereas r' scales with noise intensity with 
an exponent fairly close to 2/3. 
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FIG. 3: Effective structure-function parameters for model 
for r — 0.05 and ko = 0.6. In (a), the solid line corresponds 

0.05 + 

a(a 2 ■ 



to the linear relation a' 2 — 1.365<r 2 , and in (b) to r' 
- 2 ^ x with a ~ 0.704 and A ~ 0.6. 



The numerical results of Fig. Ufa) have been fitted with 
the analytic approximation given by Eq. (J5J) with renor- 
malized parameters, leading to a quite good agreement. 
This is shown in Fig.^Ja) for four different noise intensi- 
ties. Similarly to what occurs in the chemical model, we 
can see how noise excites a peak in the structure func- 
tion (in this case at k = ko) even within the homogeneous 
regime. This peaks increases in size as noise intensity in- 
creases, both in an absolute way and relative to the value 
of the structure function at k = 0. However, for large 
noise intensities disorder kicks in again, and the peak 
becomes less pronounced. This effect can be seen quanti- 
tatively, in terms of the dependence of AS on a 2 , for two 
different distances to threshold in Fig. Efb), where the 
numerical results are again compared to the theoretical 
prediction coming from the renormalized version of JSJ, 
which in this case gives 



AS(cr 2 ) 



r'{k§ + r') 



(7) 



The agreement between the numerical results and the 
theoretical prediction is rather satisfactory. In all cases, 
a clear enhancement of spatial coherence (in terms of spa- 
tial structures with wavenumber fc max — kg) is evident. 

In conclusion, we have demonstrated that spatiotem- 
poral noise is able to extract and enhance spatial coher- 
ence from deterministically homogeneous nonlinear me- 
dia. This effect constitutes a pure spatial analogue of 
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FIG. 4: (a) Structure function exhibited by the Swift- 
Hohenberg model Q for r = 0.005, fco = 0.6 and four differ- 
ent noise intensities, listed in the legend, (b) AS versus noise 
intensity for two different values of the control parameter r. 
Squares correspond to r = 0.005 and circles to r = 0.05. In 
the two plots, lines are the fitted function J7J with the renor- 
malized parameters given in Fig. [3] The simulations have 
been performed again in a square lattice of 64x64 cells with 
spacing Ax = 0.5. 



temporal coherence resonance, and relies on the gener- 
ation of spatial noisy precursors of pattern-forming in- 
stabilities. This phenomenon is different from what hap- 
pens in noise-induced pattern formation [28|, where spe- 
cial types of multiplicative noise induce (or displace) 
the pattern-forming transition, destabilizing the homoge- 
neous state. Here that state is still stable, and the noise 
can be simply additive. Additionally, our results lead 
us to expect that spatial analogues of both stochastic 
resonance and other mechanisms of coherence resonance 
(such as the one that takes place in excitable systems Q) 
could be found. The latter perspective, given the ubiq- 
uity of excitable systems in all areas of science is in 
our opinion very attractive. Excitable neural tissue, for 
instance, combines the features of being highly noisy and 
intrinsically spatially extended |2^| . Elucidating whether 
noise is able to enhance the spatial coherence of the sys- 
tem in this context would be extremely interesting. 

We thank L. Schimansky-Geier for useful comments. 
This research was supported by the Ministerio de Cien- 
cia y Tecnologfa (Spain) and FEDER under projects 
BFM2000-0624, BFM2001-2159, and BFM2002-04369. 
J.G.O. is partially supported by the NSF IGERT Pro- 
gram of Nonlinear Sciences (Cornell). 



* Electronic address: jordi.g.ojalvo@upc.es 
[1] K. Wiesenfeld and F. Moss, Nature 373, 33 (1995); L. 
Gammaitoni, P. Hanggi, P. Jung and F. Marchesoni, Rev. 
Mod. Phys. 70, 223 (1998); V.S. Anishchenko, V. As- 
takhov, A.B. Neiman, T. Vadivasova, and L. Schimansky- 
Geier, Nonlinear Dynamics of Chaotic and Stochastic 
Systems (Springer, Berlin, 2002). 

G. Hu, T. Ditzinger, C. Z. Ning, and H. Haken, Phys. 
Rev. Lett. 71, 807 (1993). 
[3] W.J. Rappel and S.H. Strogatz, Phys. Rev. E. 50, 3249 



[2] 



(1994). 

[4] A.S. Pikovsky and J. Kurths, Phys. Rev. Lett. 78, 775 
(1997). 

[5] B. Lindner, J. Garcia- Ojalvo, A. Neiman, and L. 

Schimansky-Geier, unpublished (2003). 
[6] L. I and J.-M. Liu, Phys. Rev. Lett. 74, 3161 (1995). 
[7] G. Giacomelli, M. Giudici, S. Balle and JR. Tredicce, 

Phys. Rev. Lett. 84, 3298 (2000). 
[8] D.E. Postnov, S.K. Han, T.G. Yim, and O.V. Sosnovt- 

seva, Phys. Rev. E 59, 3791 (1999). 
[9] D. Wilkowski, J. Ringot, D. Hennequin, and J.C. Gar- 

reau, Phys. Rev. Lett. 85, 1839 (2000). 
[10] H.G. Gu, M.H. Yang, L. Li, Z.Q. Liu, and W. Ren, Neu- 

roreport 13, 1657 (2002). 
[11] J. Garcfa-Ojalvo and J.M. Sancho, Noise in Spatially Ex- 
tended Systems (Spriger-Verlag, New- York, 1999). 
[12] S. Alonso, I. Sendiha-Nadal, V. Perez-Munuzuri, J.M. 

Sancho, and F. Sagues, Phys. Rev. Lett. 87, 078302 

(2001). 

[13] J.F. Lindner, B.K. Meadows, W.L. Ditto, M.E. Inchiosa 

and AR. Bulsara, Phys. Rev. Lett. 75, 3 (1995); H.S. 

Wio, Phys. Rev. E, 54, R3075 (1995). 
[14] S.K. Han, T.G. Yim, D.E. Postnov, and O.V. Sosnovt- 

seva, Phys. Rev. Lett. 83, 1771 (1999); A. Neiman, L. 

Schimansky-Geier, A. Cornell-Bell, and F. Moss, Phys. 

Rev. Lett. 83, 4896 (1999); C.S. Zhou, J. Kurths, and B. 

Hu, Phys. Rev. Lett. 87, 098101 (2001). 
[15] B. Lindner and L. Schimansky-Geier, Phys. Rev. E 60, 

7270 (1999). 
[16] K. Wiesenfeld, J. Stat. Phys. 19, 25 (1985). 
[17] L. Omberg, K. Dolan, A. Neiman, and F. Moss, Phys. 

Rev. E 61, 4848 (2000). 
[18] A. Neiman, P.I. Saparin and L. Stone, Phys. Rev. E 56, 

270 (1997). 

[19] I.Z. Kiss, J.L. Hudson, G.J.E. Santos, and P. Par- 
mananda, Phys. Rev. E 67, 035201 (2003). 

[20] I. Rehberg, S. Rasenat, M. de la Torre Juarez, W. Schopf, 
F. Horner, G. Ahlers, and HR. Brand, Phys. Rev. Lett. 
67, 596 (1991). 

[21] M. Wu, G. Ahlers, and D.S. Cannell, Phys. Rev. Lett. 

75, 1743 (1995). 
[22] M. Hoyuelos, P. Colet, and M. San Miguel, Phys. Rev. E 

58, 74 (1998). 

[23] IR. Epstein and J. A. Pojman, An Introduction to Non- 
linear Chemical Dynamics (Oxford University Press, New 
York, 1998). 

[24] A. Sanz-Anchelergues, A.M. Zhabotinsky, I.R. Epstein, 
and A.P. Munuzuri, Phys. Rev. E. 63, 056124 (2001). 

[25] The homogenous component of the spatial profile in this 
system gives rise to a spike of s(k) at k = 0, which has 
been extracted in the present results. The final value of 
s(0) used here is then obtained by extrapolating to k — 
a fitting function defined for k > 0. 

[26] M.C. Cross and PC. Hohenberg, Rev. Mod. Phys. 65, 
851 (1993). 

[27] J.S. Langer, Annals Phys. 65, 53 (1971). 

[28] J. Garcfa-Ojalvo, A. Hernandez-Machado and J.M. San- 
cho, Phys. Rev. Lett. 71, 1542 (1993); J.M.R. Parrondo, 
C. Van den Broeck, J. Buceta, and J. de la Rubia, Phys- 
ica A 224, 153 (1996); J. Buceta, M. Ibafies, J.M. Sancho, 
and K. Lindenberg, Phys. Rev. E 67, 021113 (2003). 

[29] D. Kleinfeld, KR. Delaney, M.S. Fee, J.A. Flores, D.W. 
Tank, and A. Gelperin, J. Neurophys. 72 1402 (1994). 



